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We study spatial structures of anomalously localized states (ALS) in tail regions at the 
critical point of the Anderson transition in the two-dimensional symplectic class. In order to 
examine tail structures of ALS, we apply the multifractal analysis only for the tail region of 
ALS and compare with the whole structure. It is found that the amplitude distribution in the 
tail region of ALS is multifractal and values of exponents characterizing multifractality are the 
same with those for typical multifractal wavefunctions in this universality class. 
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1. Introduction 

Disorder induced metal-insulator transition at abso- 
lute zero temperature, namely, the Anderson transition, 
has been extensively studied since the seminal paper by 
Anderson. 1 Although the vast amount of knowledge on 
this physical phenomenon has been accumulated for half 
a century, there still exist fundamental problems on prop- 
erties of non-interacting electrons in disordered systems. 
One of such unsettled questions is anomalously local- 
ized states (ALS) 2,3 in which most of amplitudes con- 
centrate on a narrow spatial region even in a metallic 
phase. The existence of ALS was analytically predicted, 
and the nature of ALS has been numerically and analyti- 
cally studied. 4-17 Most of previous works on ALS focused 
on the amplitude distribution function, which mainly re- 
flects peak structures of ALS. Falko et al. 5 and Uski et 
al. have predicted that the spatial correlation function 
of amplitudes has a long tail. This implies that ALS are 
not truely localized. The tail structure of ALS is crucial 
for understanding transport properties of metals. 

Recently, we have shown that ALS also exist at the 
critical point of the Anderson transition with a finite 
probability even in infinite systems. 1 It is well known 
that wavefunctions at metal-insulator transition point 
show multifractality reflecting the absence of character- 
istic lengths in critical wavefunctions. 2 1-23 Because of 
a confining length of ALS, ALS wavefunctions are not, 
however, multifractal. The finite probability of ALS im- 
plies that ALS contradict well known critical properties 
such as universality, the scaling concept, and the idea of 
renormalization. Our results suggest that these critical 
properties should be considered for typical critical states. 
It is natural to suppose that the critical level statistics 
is also influenced by ALS. Our previous works, however, 
show that ALS do not contribute to the level statistics. 20 
To understand this fact, it is quite important to reveal 
the tail structure of ALS wavefunctions. 

In this paper, we investigate spatial structures of ALS 
wavefunctions at the critical point of the Anderson tran- 
sition in the two-dimensional symplectic class. It is found 
that the amplitude distribution of ALS in their tail re- 
gion actually shows multifractality and values of expo- 



nents characterizing multifractality are the same with 
those for typical multifractal wavefunctions for this uni- 
versality class. This paper is organized as follows. In §2, 
we give a quantitative definition of ALS at the critical 
point based on the idea that ALS at criticality do not 
show multifractality. In §3, the basic multifractal analy- 
sis is remained. We briefly explain, in §4, the employed 
SU(2) model which belongs to the symplectic class and 
a numerical method to obtain eigenstates of systems de- 
scribed by this model. Numerical results and conclusions 
are given in §5. 

2. Definition of Anomalously Localized States 

In order to study ALS, we employ a definition of ALS 
proposed in ref. 18. This definition is based on the idea 
that ALS are not multifractal as a consequence of their 
localized nature. At first, we introduce the box-measure 
correlation function G q (l,L,r) defined by 24,25 



N b N br 



b r 



where fi m = £ iei)(i) |V>i| 2 and fi br(l) = £ ie6r(i) lV>i| 2 are 
box measures for wavefunction amplitudes ipi, in a box 
b(l) of size I and in a box b r (l) of size I fixed distance 
r — I away from the box 6(7), respectively. (or N br ) is 
the number of boxes b(l) [or b r {l)\, and the summation 
Tib ( or Tb ) 1S taken over all boxes b(l) [or b r (l)} in 
the system of size L. If a wavefunction is multifractal, 
G q (l,L,r) should behave as 25 



G q (l,L,r) oc i*(«)L-v(«) r -*(«), 



(2) 



where x(q), y(q), and z(q) are exponents describing mul- 
tifractal correlations of the amplitude distribution. This 
relation is sensitive to ALS as demonstrated in ref. 18 
and then suitable for defining ALS. To find the / and r 
dependences of G q (l, L, r), we concentrate on the follow- 
ing functions, 



Q q (l) = G q (l,L,r = l)ocl x ^- z ^, 



and 



R q (r)=G q (l = l,L,r) 



(3) 



(4) 
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In order to quantify non-multifractality of a specific 
wavefunction, it is convenient to introduce variances 
Var(logQ2) and Var(logi?2) from the linear functions of 
log/ and logr, log<32(0 = [x(2) — z(2)]log/ + cq and 
log_R 2 (f) = — z(2)logr + c^, respectively, calculated by 
the least-square fit. From these variances, we define a 
quantity T as 



T = AVar(logQ 2 ) + Var(logi? 2 ) 



(5) 



where A is a factor to compensate the difference between 
average values of Var(logQ2) and Var(logi? 2 )- Using T 
given by eq. (5), the quantitative and expediential defi- 
nition of ALS at criticality is presented by 



r > r*, 



(6) 



where T* is a criterial value of T to distinguish ALS from 
multifractal states. 

3. Multifractal Analysis 

In this section, we give some definitions of basic quan- 
tities and exponents used in the multifractal analysis. At 
first, we introduce a quantity Z(q) defined by 



(7) 



Here, symbols and notations in eq. (7) have the same 
meanings of those in eq. (1). For a multifractal wave- 
function, Z q (l) obeys a power law 

Z q (l) <* l T(q \ (8) 

where r(q) is called the mass exponent. From eq. (8), the 
mass exponent is given by 

log Z q (l) 



r(q) = lim • 

2^0 



log/ 



(9) 



Using r(q), the generalized dimension D q is defined as 

r(q) 



(10) 



particulary, D q for q = 2 (thus D 2 ) is called the correla- 
tion dimension. 

The multifractal spectrum /(a) is defined by the Leg- 
endre transform of r(g), i.e., f(a) = aq — r(q) 1 where the 
Lipschitz-Holdcr exponent a is defined by a — dr(q)/dq. 
The multifractal spectrum f(a) has the meaning of the 
fractal dimension of the spatial distribution of boxes 
characterized by the exponent a. The multifractal spec- 
trum f(a) and the Lipschitz-Holder a can be also calcu- 
lated directly from box measures of wavefunction ampli- 
tudes as 30 

T,b m b(l)(l) lo Sm H i)(q) 



/(«) 



lim ■ 

2^0 



log I 



and 



a = lim 



log/ 



(11) 



(12) 



respectively. Here, m b ^ (q) called the g-microscope is de- 
fined by 



m b (i)(q) 



(4(i) 



(13) 



4. System and Numerical Method 

Considering the advantage of system sizes, we fo- 
cus our attention on the Anderson transition in two- 
dimensional electron systems with strong spin-orbit in- 
teractions, in which systems have no spin-rotational sym- 
metry but have the time-reversal one. Hamiltonians de- 
scribing these systems belong to the symplectic class. 
Among several models belonging to this universality 
class, we adopt the SU(2) model because of its small scal- 
ing corrections. The Hamiltonian of the SU(2) model 26 
is given by 



(14) 



where cj (cj) is the creation (annihilation) operator act- 
ing on a quaternion state vector, R^ is the quaternion- 
real hopping matrix element between the sites i and j, 
and Si denotes the on-site random potential distributed 
uniformly in the interval [—W/2,W/2]. (Bold symbols 
represent quaternion-real quantities.) The matrix ele- 
ment is given by 

Rij — cos OLij cos PijT + sin 7^ sin fyjT 

— cos -fij sin PijT 2 + sin cos flijT 3 , (15) 

for the nearest neighbor sites i and j, and R^ = for 
otherwise. Here, r M (/1 = 0,1,2,3) is the primitive ele- 
ment of quaternions. 27 Random quantities and 7^ 
are distributed uniformly in the range of [0, 2n), and 
Pij is distributed according to the probability density 
P(/3)d/3 = sin(2/3)d/3 for < (3 < tt/2. Randomly distrit- 
buted hopping matrix elements shorten the spin relax- 
ation length which is a dominant irrelevant length scale. 
Thus, scaling corrections become very small in the SU(2) 
model. It is known that the localization length exponent 
v of this model is 2.73 ±0.02 and the critical disorder W c 
is 5.952U at E = 1.0U. 26 

Critical wavefunctions of the SU(2) model have 
been calculated by using the forced oscillator method 
(FOM) 28 extended to eigenvalue problems of quaternion- 
real matrices. Of course, the Hamiltonian eq. (14) can 
be represented by complex numbers, and we can use the 
usual FOM for complex Hermitian matrices to solve the 
eigenvalue problem. The modified FOM for quaternion- 
real matrices, however, enables us to calculate eigenval- 
ues and eigenvectors within about a half of CPU time. 29 
It should be remarked that the obtained eigenvector 
is a quaternion-real vector. This vector represents two 
physical states simultaneously, which correspond to the 
Kramers doublet. Since the amplitude distribution of 
these degenerate states are the same, we analyze one of 
the calculated Kramers doublet. 

5. Results and Conclusions 

In this section, we show the results of the multifractal 
analysis for amplitude distribution of ALS in the tail re- 
gion. At first, we look for ALS according to our definition 
of eq. (6). In this work, the criterial value is T* = 0.03. 
For this purpose, we calculate 10 4 critical wavefunctions 
for the SU(2) model at E = 1.0V and W = 5.952V by 
means of the FOM. Each eigenstate is obtained for a 
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Fig. 2. Z2U) for -0ALS (open circles) and i/'tail (filled circles). The 
dashed line shows the least square fit for Z2(l) for Vtaii- 



Fig. 1. Squared amplitudes of a wavefunction V'ALS f° r th e SU(2) 
model of the system size L = 120 (lower figure). The value char- 
acterizing non-multifractality defined by cq. (5) is T = 0.10. The 
square enclosed by solid lines represents a tail region. The upper 
figure shows squared amplitudes of the wavefunction Vtaii of this 
tail region (the subsystem size is L = 60). 

single disorder realization. Periodic boundary conditions 
arc imposed in the x and y directions in the system of 
size L = 120. 

Figure 1 shows squared amplitudes of a ALS wave- 
function V'als- The value of T of this wavefunction is 
r = 0.10. In this wavefunction, amplitudes concentrate 
in a narrow spatial region. The wavefunction appears to 
be localized for this feature. In order to investigate the 
tail structure of this wavefunction, we extract a part of 
the wavefunction within a 60 x 60 subsystem where the 
cut region is depicted in Fig. 1 by solid lines. The upper 
figure of Fig. 1 shows the extracted tail region wave- 
function V'taii- The amplitude distribution of V'taii is very 
complicated, and seems to be multifractal. To clarify this 
point, we examine below multifractality of V'taii- 

At first, we calculate Z 2 {1) given by eq. (7) for V'als 
and V'taii- The wavefunction V'taii is normalized in the 
subsystem to perform the multifractal analysis. We see 
that Z 2 (l) for V'als does not exhibit the power law as 
shown in Fig. 2 (open circles). If ALS are truly localized 
such as V'als M °c exp(— r/£), where £ is a localization 
length, Z 2 (l) for V'als should take a constant value for 
I £. However, Z 2 {1) for ALS is not constant even for 
large I as shown in Fig. 2. On the contrary, the quan- 
tity Z 2 (l) for V'taii clearly obeys a power law as shown 
by filled circles in Fig. 2. The value of the correlation 
dimension D 2 = t(2) = 1.69 ± 0.01 is very close to the 
value reported so far for typical multifractal wavefunc- 
tions belonging to this universality class. 18 This implies 
that tail regions of ALS exhibit the same multifractality 
with that of typical critical wavefunctions. The fact that 
the behavior of Z 2 {1) for V'als for Z S> 20 is similar to 
that for Vtaii shows that Z 2 (l) at large I values is domi- 
nated by the tail structure of the ALS wavefunction. 

Next, we calculate the multifractal spectra for both 



wavefunctions V'als an d V'taii- Results are shown in 
Fig. 3. Remarkably, f(a) for Vals takes negative val- 
ues for small a, which is not reasonable because f(a) is 
the fractal dimension of the distribution of boxes char- 
acterized by a. This supports that Vals does not have 
a multifractal structure. On the other hand, f(a) for 
V'taii satisfies several conditions for multifractal spectra. 
The value of ceo giving the maximum of f(a) for V'taii is 
cxq = 2.15 ± 0.02, and the minimum and maximum val- 
ues of a are estimated as a m i n = 1.22 and a max = 3.26, 
respectively. These values are very close to those for typi- 
cal multifractal wavefunctions. In fact, we calculate f{a) 
for a typical critical wavefunction with T — 0.001 which 
is the smallest value of T in all calculated wavefunc- 
tions. Obtained values are a = 2.17, a m i n = 1-22, and 
ttmax = 3.31. The characteristic values of a for V'taii are 
also close to averaged values for an ensemble of typical 
critical wavefunctions. 19 The profile of f(a) for Vtaii in 
the vicinity of a = a can be well approximated by the 
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Fig. 3. Multifractal spectra /(») for V'als (open circles) and -i/'tail 
(filled circles). The dashed line shows the parabolic approxima- 
tion eq. (16) with cxq = 2.15. 
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Fig. 4. The generalized dimension as a function of the moment q 
for i/'ALS (open circles) and Vtail (filled circles). The dashed line 
represents the parabolic approximation cq. (17) with «o = 2.15. 



parabolic form 



/(«) 



(q - af 
4{a - 2) ' 



(16) 



which is shown by dashed line in Fig. 3. These results 
ensure that Vtaii has the same multifractal amplitude 
distribution with that for typical critical wavefunctions. 
It seems that /(a) for Vals is also fitted by eq. (16). 
This possibility is, however, denied when examining the 
generalized dimension D q . 

Figure 4 shows the generalized dimension D q calcu- 
lated by the box-counting method eq. (7) for both wave- 
functions V'als and V'taii- The values of D q for two wave- 
functions largely deviate each other for \q\ 3> 1. Since 
D-oo and are equal to a max and a m - m , respectively, 
such large deviations are consistent with Fig. 3. The 
parabolic approximation of D q corresponding to eq. (16) 
is given by 



D q = 2 - q(a - 2), 



(17) 



which would be valid in the vicinity of q = 0. The gen- 
eralized dimension D q for iptnii is well approximated by 
eq. (17) with ao = 2.15 as shown by dashed line in Fig. 4. 
On the contrary, D q for f/'ALS does not follow a straight 
line near q = 0. This implies that the parabolic approx- 
imation eq. (16) does not reproduce f(a) for t/als even 
in the vicinity of a = ao- 

In summary, we investigate the spatial structure of 
ALS at the critical point of the Anderson transition. It is 
found that ALS are not truly localized and the tail struc- 
ture of ALS shows the same multifractality with that of 
typical wavefunctions. These results form the foundation 
of understanding the distribution of physical quantities 
at criticality. 
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